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O . Abstract 



The decoupling properties of the Higgs sector in the Minimal Supersymmetric Standard Model (MSSM) 
imply that a light CP-even Higgs boson discovered at the Tevatron or LHC may closely resemble the 



^ • 

Oh 
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Standard Model (SM) Higgs boson. In this paper, we investigate how precision measurements of Higgs 



properties at a Linear Collider (LC) can distinguish between a CP-even Higgs boson of the MSSM and 
the SM Higgs boson. We review the expected theoretical behavior of the partial widths and branching 



H " ratios for decays of the neutral MSSM Higgs bosons with significant couplings to the W and Z bosons, 
including the leading radiative corrections to the mixing angle a and tan /3-enhanced vertex corrections. 
The general expectation is that the Higgs couplings to W~^W~ , ZZ, cc and ti should quickly approach 
their SM values for increasing CP-odd Higgs mass ruA, while the couplings to bb and t^t~ do so more 
slowly. Using the expected experimental and theoretical accuracy in determining SM branching ratios 
and partial widths, we demonstrate the sensitivity of measurements at the LC to variations in the MSSM 
parameters, with particular attention to the decoupling limit. For a wide range of MSSM parameters, the 
LC is sensitive to niA ~ 600 GeV almost independently of tan (3. For large values of tan (3 and some specihc 
choices of MSSM parameters [e.g., Atfj, < and \At\ — \fj,\ — 0{Ms)], one of the CP-even Higgs bosons 
can be SM-like independent of the value of tha- In the case of large deviations from the SM, we present 
a procedure using Higgs coupling measurements to extract the supersymmetric correction to the relation 
between the b quark mass and Yukawa coupling. 

'Electronic addresses: carena@fnal.gov, liaber@scipp.ucsc.edu, logan@fnal.gov, mrenna@physics.ucdavis.edu 



1 Introduction 



The radiative corrections to Higgs boson masses and couplings in the minimal super symmetric extension 
of the Standard Model (MSSM) have been investigated thoroughly using different theoretical approaches. 
Derived quantities such as Higgs boson production cross sections, partial widths and branching ratios 
(BRs) are predicted to a high level of precision for any given set of MSSM parameters. While at least one 
of the neutral Higgs bosons of the MSSM has a coupling to W and Z bosons similar in magnitude to a 
Standard Model (SM) Higgs boson, some of its properties can differ from those of the SM Higgs boson of 
the same mass. Nevertheless, over a significant region of parameter space, the deviation of the couplings of 
this SM-lik^ Higgs boson from the corresponding couplings of the SM Higgs boson is small and approaches 
zero in the so-called "decoupling limit" of the model ^ . 

Experiments at the Fermilab Tevatron Q and CERN LHC |^ will be sensitive to the Higgs bosons 
of the SM and MSSM. The Tevatron can discover a SM-like Higgs boson for most choices of MSSM 
parameters if enough data can be accumulated, the detectors perform as expected, and systematic errors 
are demonstrably small. The LHC can discover at least one MSSM Higgs boson over all of the MSSM 
parameter space and several Higgs bosons are likely to be discovered in a significant region of the 
parameter space. For moderate values of the ratio, tan/?, of the two Higgs vacuum expectation values,0in 
the range 2.4 < tan/3 < 8 [ 2.4 < tan/3 < 17] for a CP-odd Higgs mass of niA = 250 GeV [ruA = 500 GeV], 
only one SM-like Higgs boson will be visible at the LHC ||3|,|5). If more than one Higgs boson is observed 
at the hadron colliders, then additional precision measurements of Higgs decay properties can determine if 
the Higgs bosons originate from a two-Higgs doublet sector of a supersymmetric (SUSY) model (as in the 
MSSM) or from a different model. If only one Higgs boson is observed, such precision measurements could 
indicate whether any additional Higgs structure exists. Of course, the MSSM also contains supersymmetric 
particles, which can be discovered at the Tevatron and/or the LHC. However, these particles can be heavy, 
so that only a part of the spectrum may ultimately be observable at the LHC, or the interpretation of the 
data as SUSY particle production may be ambiguous. Hence, precision measurements of the properties 
of the Higgs sector can provide crucial supporting evidence for the MSSM. In this study, we explore the 
potential of a future e~^e~ linear collider (LC) to explore the MSSM Higgs sector in regions of MSSM 
parameter space with very different behaviors. 

If a Higgs boson couples to Z bosons with SM-like strength, then its mass can be determined to high 
precision at the LC through measurements of the recoil mass spectrum against a Z boson. With 500 fb~^ 
^In this paper, a SM-like Higgs boson always refers to the neutral Higgs boson of the MSSM with gfi^.yy ^ ^-^ghsMVV j 
where V = W oi Z. The decay properties of this Higgs boson may be quite different than those of a SM Higgs boson in the 

parameter regime away from the decoupling limit. 

^The absence of a Higgs boson discovery at LEP implies that the range 0.5 < tan/9 < 2.4 is excluded at 95% confidence 

level [§. 
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of data at y/s = 350 GeV, a precision of 40-90 MeV can be achieved for Higgs masses between 120 and 
180 GeV [^,^. With the mass of the Higgs boson so constrained, the theoretical predictions for many of 
the Higgs branching ratios and partial widths (i.e., the Higgs couplings) will be known to great accuracy. 
As a result, precision experimental measurements of Higgs branching ratios and partial widths may allow 
one to discriminate between the SM and new physics of electroweak symmetry breaking. A number of 
recent studies have evaluated how precisely the branching ratios and couplings of the SM Higgs boson 
can be measured at the LC |7-13|. In this paper, we exploit the results of these studies to examine the 



potential of Higgs boson BR and partial width measurements at the LC to distinguish the SM-like Higgs 
boson of the MSSM from the SM Higgs boson. We focus on choices of MSSM parameters that exhibit a 
significant variation in the approach to the decoupling limit. This allows us to identify parameter regions 
where hadron collider measurements will most likely not be able to distinguish a MSSM from the SM Higgs 
boson, and LC measurements can be essential for this purpose. In addition, one new observation is that 
there are regions where decoupling occurs at fairly low values of so that precision measurements will not 
be able to reveal much about the MSSM parameter space structure. In such cases, direct measurements of 
the properties of the heavier Higgs bosons will be necessary to elucidate the true nature of the Higgs sector. 
Even if the existence of a MSSM Higgs boson is established, it may be challenging to extract the underlying 



MSSM parameters [14|. Nevertheless, in some regions of MSSM parameter space, we demonstrate that 
the SUSY vertex corrections to the Higgs boson couplings to bottom quark and tau pairs can be extracted 
from branching ratio and Higgs coupling measurements]^ 

This paper is organized as follows. In Sec. |2|, we outline the features of the MSSM Higgs sector relevant 
for our analysis. In particular, we review and expand on the details of the radiative corrections that 
contribute significantly to the Higgs couplings. In Sec. ^, we analyze the behavior of the Higgs partial 
decay widths in various regions of MSSM parameter space. We consider three benchmark scenarios that 
lead to very different behaviors in the MSSM Higgs sector. In Sec. ^, we combine the expected experimental 
capabilities of the LC to measure Higgs boson BRs and couplings with our theoretical results. For the 
benchmark scenarios, we quantify the ability of the LC to distinguish the lightest MSSM Higgs boson from 
the SM Higgs boson. In Sec. |5|, we describe how the SUSY Yukawa correction to the Higgs boson couplings 
to bottom quark pairs, Af,, can be extracted from Higgs measurements. In Sec. ^ and ^ we also evaluate 
the impact of the theoretical uncertainties on our analysis. Finally, Sec. ^ contains our conclusions. Some 



preliminary results of this work were presented in Ref. |17]. 



^The distinctive effects of SUSY vertex corrections to Higgs-fermion pair couplings on brandling ratio measurements have 
also been explored in Refs. [ p^[l6| . 
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2 The MSSM Higgs sector 



In this section, we review those properties of the MSSM Higgs sector relevant to our analysis. At tree 
level, the masses and couplings of the MSSM Higgs bosons are determined by two parameters, which 
are conveniently chosen to be the mass of the CP-odd Higgs boson, tha, and the ratio of the vacuum 
expectation values of the two neutral Higgs fields, tan/3. Radiative corrections to the MSSM Higgs sector 
introduce significant dependence on other MSSM parameters (for a review see Ref. Q). These radiative 
corrections have been analyzed extensively in the literature P,[I^-p5[. 

The two main sources of radiative corrections to the couplings of the MSSM Higgs bosons are: (i) the 
radiative corrections to the Higgs squared-mass matrix p"8|-^ , which give rise to corrections to an effective 
CP-even Higgs mixing angle a, and (ii) vertex corrections to the Higgs- fermion Yukawa couplings |27,^9|-45|. 
In this paper we examine the effects of these two types of corrections in the CP-conserving MSSM.0 

The squared-mass matrix for the CP-even neutral MSSM Higgs bosons h and H (where nih < mn) is 
given by: 

j^2^(^n -^12^^/ + -("^a + "^I)s/3C/3 \ ^_^2 (2 1) 

\Ml2 MI2 j \-{m\ + ml)spCp rn\cl + mlsl j 

where SM"^ is a consequence of the radiative corrections. At tree level, one obtains rrih < mz\ cos 2/3| < mz- 
Such a light h is essentially ruled out by searches at LEP2j^ However, once radiative corrections to the 
squared-mass matrix are included, the theoretical upper bound on ruh is raised substantially. For a fixed 
value of tan/3 and a specified set of MSSM parameters, rrih grows with increasing niA and reaches an 
asymptotic value m™'^^(tan/3) in the limit of large rriA- If tan/3 is now allowed to vary (while holding all 
other free parameters fixed), m™'^^(tan/3) increases with tan/3 and typically^ reaches an asymptotic value 
j^max £qj, tan/3 > 10. For large values of tan/3, ruh ~ m™'^-^ and mn ~ ttia for niA > m^^^. Conversely, if 
niA < m^^^ then rrih ~ rriA and mn ~ m™^^. 

At the LC, a light SM-like Higgs boson will be produced through Higgsstrahlung [e'^e~ Z* ^ Zh], 
WW fusion [e+e~ W*W*iyi> huv] and ZZ fusion [e+e~ Z*Z*e^e~ /le+e^]. The cross sections 
for all the above processes depend on the couplings of the Higgs boson to vector boson pairs. In the MSSM, 
the couplings of h [H] to vector boson pairs are given by sin(/3 — a) [cos(/3 — a)] times the corresponding 
SM Higgs coupling. The decoupling limit corresponds to tua ^ rnz, in which case sin(/3 — a) ~ 1, and the 
properties of h approach those of the SM Higgs boson. The CP-even Higgs squared-masses obey the sum 
*The MSSM Higgs sector automatically conserves CP at tree-level, although non-trivial CP-violating effects can enter at 
one-loop (due to complex MSSM parameters) and be phenomenologically significant |^,^. In this paper, we assume that 

the one-loop CP-violating effects are absent. The CP-violating case will be addressed elsewhere. 
'^The current MSSM Higgs mass limits are > 91.0 GeV and ruA > 91.9 GeV 

®ln some regions of MSSM parameter space at large tan fJ, radiative corrections to the Higgs-bottom quark Yukawa coupling 
can yield large negative loop corrections to m™^''(tan/3), so that the latter begins to decrease for tan/3 > 10. 
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rule 11,11 

mjf cos'^ {/3 - a) + ml sin^ip- a) = [m'^''''{tan p)f . (2.2) 

In particular, combining this sum rule with the large tan /3 behavior of the Higgs masses noted above 
implies that H is the SM-like Higgs boson for mA < mf^^^ and large values of tan (3, while h is the SM-like 
Higgs boson for mA > mf^^^. Note that the decoupling limit implies that the latter holds for any value of 
tan (3. Most of the analysis of this paper will focus on the case where h is the SM-like Higgs boson and 
mA > mf^^^. 

Along with tan /3, the CP-even Higgs mixing angle a determines the Higgs boson couplings to fermions. 
In particular, relative to their SM values, the couplings of h [H] to down-type fermions are multiplied 
by — sin a/ cos /3 [cos a/ cos /?], and those of h [H] to up- type fermions are multiplied by cos a/ sin/? 
[sin a/ sin /3] . Thus radiative corrections to a can have significant effects on the Higgs boson couplings 



to fermions. The mixing angle a which diagonalizes the mass matrix in Eq. 2.1 can be expressed as: 



M 



12 



22 



(2.3) 



where Sa = sin a and Cq = cos a. Note that if 7Wx2 ~^ 0' then either sin a (if A^fi > -^22) 
(if Mil < -^22)- At tree level (see Eq. |2l|), MI2 

is small for small and/or large tan/3, but 



cos a 



it cannot vanish. This is no longer true after including radiative corrections, which can be of the same 
order as the tree level value for small values of mA and large tan (3. In particular, the radiatively-corrected 
value of M12 exhibits a widely varying behavior as a function of the MSSM parameters. The radiative 
corrections to A^^, including dominant corrections coming from the one- loop top and bottom quark and 
top and bottom squark contributions plus the two- loop leading logarithmic contributions, are given to 



Oiht,ht) by 



6Mii ~ -fi xt 



2 2 ^ 



32^"^ 



1 + cii In ^ 



mt 



2„2%^2 



1 + ci2 In ^ 
\ mi 



^22 



sp In 



m? 



5M 



12 



87r2 

+^^SpXtat{l2 - xtat) 

- ^t^^R ^ 2 

-/^a^t^^(6 - xtatjsp 



1 + ^21 In -4 



Mf 



mt 



1 + C21 In ^ 
\ mi , 



1 + C31 In ^ 
\ mi 



-ihfv^ 2 



-s ^h^^ 2 



1 + C22 In ^ 



1 + C32 In ^ 




where = sin/3, cp = cos/3, and the coefficients Cjj are: 



32<7i 



with {tu,ti2,t2i,t22,t3i,t32) = (12,-4,6,-10,9,-7) and (611,612,621,^22,^31,^, 



(2.4) 



(2.5) 



32) 



:-4, 12, 2, 18, -1,15). 

Above, ht and hi, are the top and bottom quark Yukawa couplings [see Eqs. |2.6| - |2.8| , gs is the strong QCD 
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coupling, V = 246 GeV is the SM Higgs vacuum expectation value, and M| = ^{M^^ + M?^) is the average 
squared top squark mass.[| The 6A4fj also depend on the MSSM parameters At, A^, and /i that enter 
the off-diagonal top-squark and bottom-squark squared-mass matrices. We employ the following notation: 
/i = iijMs, at = At/Ms, ttb = Aij/Ms and xt = Xt/Ms, where Xt = At — /icot/3. Note that the leading 
radiative corrections to A^i2 depend strongly on the sign of fj,Xt and the magnitude of At. For the scenarios 
we shall consider, with <^ Q, the combination Atfi < [Atfi > 0] can lead to a suppression [enhancement] 
of Aii2 = —{fn\ + fn'z)spcp + 5M.'i2^ ^"^^ hence, to a suppression [enhancement] of the coupling of the 
SM-like Higgs boson to b quarks and r leptons. 

Our numerical calculation of the radiative corrections to the Higgs masses and mixing angle is based on 
the results of Ref. |36| and incorporates the renormalization group improved one- loop effective potential. 



plus the non-logarithmic two-loop contributions of the Yukawa vertex corrections for top and bottom 
quarks. The Yukawa vertex corrections modify the effective Lagrangian that describes the coupling of the 
Higgs bosons to the third generation quarks:^ 



- C^s = [{h + 6h)bRWaQi + {ht + 6ht)iRQlHi\ + AhttRQim* + AhbRQlH^* + h.c. , (2.6) 
resulting in a modification of the tree-level relation between ht [hf,] and rrit [mi,] as follows: 

r., = cos /? (^1 + — + j = ^ cos /5(1 + A,) , (2.7) 

mt = -^s.np[l + — + j = ^ sm /3(1 + A,) . (2.8) 

The dominant contributions to A;, are tan /3-enhanced, with A^ ~ (A/ib/Zi;,) tan /3; for tan/3 ^ 1, Sh^/hb 
provides a small correction to A;,. In the same limit, A^ ~ 6ht/ht, with the additional contribution of 
(Aht/ht) cot /3 providing a small correction.^ Explicitly, one finds that for tan/3 > 1 |23,|4|,||] 



A,, ~ 



^/.M, I (Ml ,Ml,Ml) + r^^^At ml ,Ml,„^) 



tan 13 , (2.9) 



- -^AtM,IiMl,Ml,M^) - ^/x^/(M^^,M^^,^^), (2.10) 

where Og = g^/^ir, Mg is the gluino mass, Mj^^^ are the bottom squark masses, and smaller electroweak 
corrections have been ignored. The loop integral /(a^, b'^,<?) is of order l/max(a^, b"^ ,(?) when at least one 
of its arguments is large compared to m^; for the explicit expression see Ref. [^. These one- loop Yukawa 



'^Eq. is derived under the assumption that |M? - M? |/(M?^ ) 1- The approximate forms of Eq. are sufficient 

to provide insight on the dependence of the radiatively-corrected Higgs masses and couplings on the MSSM parameters. The 

numerical work of this paper employs more exact expressions for the Ifiggs squared-mass matrix elements as noted below. 
*The Higgs couplings to leptons and to first and second generation quarks can be treated similarly. 

^Because the one-loop corrections 5hh, Ahb, Sht and A/it depend only on Yukawa and gauge couplings and SUSY particle 
masses, they contain no hidden tan/3 enhancements |50[. 
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vertex corrections enter indirectly as two-loop effects in the Higgs squared-mass matrix elements via the 
dependence on ht and /i^. 

We have noted earlier that in the decoupling limit sin(/3 — a) = 1 [or equivalently cos(/3 — a) = 0], in 
which case the couplings of h are identical to those of the SM Higgs boson. This limit is achieved when 
m-A ^ mz- This behavior, which is easy to verify for the tree- level expressions, continues to hold when 
radiative corrections are included. However, the onset of decoupling can be significantly affected by the 



radiative corrections, as we now discuss. From Eq. 2.3, one easily obtains: 
cos(/3 — q) 



{Mil - -^22) sin 2p - 2MI2 cos 2/3 



2{mjj — ml) sin(/3 



a 



m| sin 4/? + {6Mli - SMI2) sin 2/? - 25MI2 cos 2/3 



a 



Since 6Mfi ~ 



0{m\), and 



ml 



cos(/3 — a) 



2(m^ — m\) sin(/3 
m\ + 0{m\), one obtains for mz 
m\ sin 4/3 



(2.11) 



2m\ 



+ 



mr 



m^, 



where 



1 + 



6M^ 



22 



5M^ 



12 



2mi cos 2/3 



my sin 2/3 



(2.12) 



(2.13) 



Eq. 2.12 exhibits the expected decoupling behavior for m^ ^ mz- However, Eq. 2.11 exhibits another way 
in which cos(/3 — a) = can be achieved — simply choose the super symmetric parameters (that govern the 



Higgs mass radiative corrections) such that the numerator of Eq. 2.11 vanishes. That is 



2m| sin 2/3 = 2 6Mj2 - tan 2/3 (^6M 



6MI2 



(2.14) 



sin 2/3 ~ —2/ tan/3, one can determine the 



Note that Eq. 2.14 is independent of the value of m^. For a typical choice of MSSM parameters, Eq. 2.14 
yields a solution at large tan /3 — by approximating tan 2/3 : 
value of /3 at which the decoupling occurs: 

5Mli + 5MI2 



tan (3 



2m| 



5MI2 



(2.15) 



The explicit expressions for 6A4fj quoted in Eq. 2.4 confirm that the assumption of tan /3 ^ 1 used to 



derive this result is a consistent approximation because 6Aii2 is typically small. We conclude that for the 
value of tan /3 specified in Eq. |2.15 , cos(/3 — a) = independently of the value of mA- We shall refer to this 



phenomenon as m^-independent decoupling. From Eq. 2.4, it follows that explicit solutions to Eq. 2.14 



depend on ratios of SUSY parameters and so are insensitive to the overall SUSY mass scale, modulo a 
mild logarithmic dependence on Ms/mt. 



Eq. [2.14 is equivalent to the condition c = [see Eqs. 2.12 and ^.1 



7 



The introduction of the radiatively-corrected value for the CP-even Higgs mixing angle a affects the 
MSSM Higgs boson couplings to all down-type fermions (and likewise to all up-type fermions) in the same 
way. In particular ,[3 BR(6)/BR(r) = r(6)/r(r) oc gl^Jgl^^ cx ml/ml |l5|,|l|,|3|,|5l| . However, the Yukawa 
vertex corrections enter directly in the couplings of fermions to the Higgs bosons. These corrections can 
be understood as a modification of the relation between the fermion Yukawa coupling and its mass, as 



exhibited in Eqs. 2.7- 2. J. After including the dominant corrections, the CP-even Higgs boson couplings to 



h and r are modified relative to the SM coupling, ghsmff ~ 9^f /'^''^w , as follows [43|: 



1 



9hbb 



9Hbb 



gtUb sm a 
2mw cos (3 1 + Ab 

grrib cos a 1 



Sh 

1 — A5 cot a cot (3 H (1 + cot a cot f3) 

h 

Sh 

1 + Af, tan a cot /3 + — — (1 — tan a cot /3) 
hb 



2mvycos/3 1 + Af, [ ' " ^ ' hb ^ ""'J ^ ^ 

At large tan /3, terms involving A5 oc tan f3 [Eq. |2.1C1(| provide the dominant corrections to the neutral Higgs 



couplings to hb. The corrections proportional to 5hb/hb [see Eqs. 2^-2^ and the discussion that follows] 
are never tan /3-enhanced and are therefore numerically unimportant. The r couplings are obtained from 
Eq. 2.16| by replacing rub, Ab and Shb with m,-. A,- and 6hr, respectively. Note that A,- and 5h-r arise in 
analogy with A;, and 6hb from SUSY particle loops involving the leptonic sector. In particular, at large 



tan/3 [16,27,4C,4] 



A^ ~ 



^M^^^IiMl,Ml,Ml 



tan (3 , 



(2.17) 



where 02 = 5^/47r and ai = g'^/47r are the electroweak gauge couplings. Since corrections to hr are 
proportional to qi and 02, we expect |At-| <^ \ Ab\. Nevertheless, we shall formally keep the A^ corrections 
in our analysis, although they will have negligible effect in our numerical results. 

To see how the decoupling limit is achieved for the hbb (and hrr) couplings, note that we can write: 



sma 
cos (3 



sm{(3 — a) — tan /3 cos(/3 — a) . 



Working to first order in cos(/3 — a), and usingP^ 

tan a tan /3 = — 1 + (tan/3 + cot /3) cos(/3 — a) + O ^cos^(/3 

it follows that 



a 



ghbb — 9hsMbb 



1 — (tan /3 + cot /3) cos(/3 — a) ( sin^ /3 



Ab - Shb/ hi 
1 + Aft 



(2.18) 



(2.19) 



(2.20) 



Two points are particularly noteworthy. First, for ^ rnz, decoupling is achieved since cos(/3 — q) ~ 
0{m'^z/m\) [Eq. ^.12 |. However, because A;, oc tan/3, decoupling can be delayed until m\ ^> m| tan/3 [44|. 

^^We use the notation BR(6) = BR(/!, — » bb) and r{b) = r{h bb), and similarly for other Higgs decay final states. 
^^In the decoupling limit {niA S> mz), Eq. 2.1i implies that that |(tan P + cot P) cos(/3 — q)| <C 1 for all values of tan p. 
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Second, as noted above, if Eq. 2.14 is satisfied (or equivalently if c = [Eq. p.l3| ), decoupling is achieved 
independently of the value of mA_- One also obtains a similar expression for g^rT by replacing Af, and 
with At- and 6hr in Eq. 2.20| . Since |At-| <^ \^b\, it follows that the SM expectation, g\i,i,/ gf^T-j. oc m'^/m'^ 
is violated except in the decoupling limit [15, 16, 4^,^. 

Another limiting case of interest is one where sin(/3 — a) is close to zero. This limit can be reached for 
values of itia ^ m^^^ and large tan (3. In this limit, the H couplings to the W and Z bosons approach 
their Standard Model values and become relevant for our analysis. However, this is not a decoupling limit, 
and the H couplings to fermion pairs can deviate from the corresponding Standard Model couplings. This 
can be confirmed by observing that 



cos a 



cos j3 



cos(/3 — q) + tan (3 sin(/3 



a) 



(2.21) 



But, for large tan/?, it is possible to have tan/3sin(/3 — a) ~ 0{\) even in the limit of small sin(/3 — a) 
Working to first order in sin(/3 — a) and using 



cot a tan /? = 1 + (tan /? + cot /3) sin(/3 — a) + O [s\v?{f3 — a 



it follows that 



gmb — 9?tRMbb 



1 - (tan /3 + cot (3) sin(/3 



a 



cos^ [3 



1 + 6hb/hb 



(2.22) 



(2.23) 



Note that if |(tan /? + cot /?) sin(/3 — a)| <C 1, then the Hbb coupling approaches the Standard Model value, 
even when the Yukawa vertex corrections are included. 

We next consider the CP-even Higgs boson couplings to top quark pairs. The analysis is similar to the 
one given above, and one obtains 



ghtt 



gmt 



gmt cos a 
2mw sin /? 

gmt sin a 
2mw sin (3 



1 



1 



1 A^f 

1 + At ht 

1 Aht 



(cot P + tan a] 
(cot P — cot a) 



(2.24) 



1 + At ht 

Here, it is more convenient to express our results in terms of A^ and Aht/ht, since A^ ~ 5ht/ht, while the 
corresponding contribution of Aht/ht is tan /3 suppressed [Eq. |2.8|. The Higgs couplings to charm quark 



pairs are obtained from Eq. 2.24 by replacing mt, At and Aht with rric, Ac and Ahc, respectively. Using 
cot P + tan a ~ cos(/3 — a)/ sin^ /? [see Eq. 2.19| ] , it follows that the SUSY vertex corrections to ghtt and ghcc 
are suppressed in the decoupling limit (with no enhancement in the limit of large tan /3) , and so ghtt and 
ghcc approach their Standard Model values. In the opposite limit in which sin(/3 — a) is close to zero, we 
use cot P — cot a ~ — sm{P — a)/ sin^ P [see Eq. 2.22 1 to conclude that the SUSY vertex corrections to gntt 
and gHcc are suppressed. Thus, in the numerical results presented in this paper, the SUSY Yukawa vertex 
corrections to both h and H couplings to up-type quark pairs have a negligible effect in the parameter 
regions of interest and can be neglected. 
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By considering certain ratios of Higgs-fermion couplings, one can begin to isolate various combinations 
of SUSY Yukawa vertex corrections. We introduce the notation 5^// = Qcpf f / ghsMf f = ^' -^1 ^he 
Higgs-fermion couplings normalized to their Standard Model values. From Eqs. |2.1(: and 2.24| , we obtain 



Ab-Ar 6hb _^ f 1 + Ab\ 6hr 



ghbb — 9hTT _ gmb - gHrr ^ 1 + At- hb \l + Ar J K ^2 25) 

ghtt-ghbb gmt -ghbb ~ ^ / l + AA Aht 6hb 

U + aJ ht ^ hb 

Note that dependence on the CP-even Higgs mixing angle has conveniently canceled out. At large tan f3, 
one can to first approximation keep only those terms that are tan /3-enhanced at one loop. We then obtain: 

9hbb — QhTT gubb — guTT Ab — Aj. 

— = — — 7 — — Ab , [2.2b) 

Qhtt — 9hbb gmt — ghbb 1 + A^ 

where the last step follows if | A,- 1 <^ 1 (we have already noted that | A,- 1 <^ | A;, | ) . Eq. 2.26| and the analogous 



result in which the Higgs couplings to top quarks are replaced by the corresponding couplings to charm 
quarks will be used in Sec. ^ when we discuss the extraction of A^ from Higgs coupling measurements. 

3 Behavior of Higgs decay observables 

In this section, we examine in detail the behavior of the MSSM Higgs boson partial widths. In order 
to present quantitative results, we consider three "benchmark" scenarios for the MSSM parameters that 
lead to very different behaviors of the SM-like Higgs boson of the MSSM. Our three benchmark scenarios, 
summarized in Table |l|, correspond approximately to those discussed in Ref. [31|. All MSSM parameters 
are specified at the electroweak scale. The three benchmark scenarios have the following properties: 

No-mixing scenario: The top squark mixing angle % is zero. This scenario yields the lowest value of 
m™'*^(tan/?) for given values of tan/3 and Ms- For simplicity, we define the scenarios in terms of 
-^SUSY = Mq = = Mj=), where the latter are third generation squark mass parameters. For 
^SUSY ^ frit, as is true in the scenarios considered here, MgusY — Ms [where M| = ^{M~^ + M£)]. 
Here we have chosen a large value for MsusY = 1-5 TeV in order to obtain a sufficiently large value 
of m™^^(tan/3), comparable to that obtained in the other two scenarios (the case of AfsusY = 1 TeV 
is at the edge of the region excluded by LEP2). 

Maximal-mixing scenario: The top squark mixing is chosen to give the maximal value of m™^^(tan/3) 
for given values of tan (3 and Ms- 

Large fj, and At scenario: Large radiative corrections occur to both a and A5. In particular, Aii2 can 
exhibit extreme variations in magnitude depending on the sign of Atfj, and the magnitude of Af. 
The two possible sign combinations for At and fi (for a fixed sign of Atfi) yield small differences in 
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Mass paramet 


3rs [TeV 








Benchmark 




Xt = At — fi cot f3 


-^SUSY 


Mg 


m^'^^ [GeV] 


No-Mixing 


-0.2 





At 


1.5 


1 


118 


Maximal- Mixing 


-0.2 


Vg 


At 


1 


1 


129 


Large /i and At 


±1.2 


=Fl.2(l + cot /?) 





1 


0.5 


119 



Table 1: MSSM parameters for our benchmark scenarios, and the derived maximal mass for the SM-like 
Higgs boson. 

Mi2 through the dependence of ht and hh on and Af,, respectively. The vertex correction A^ 
is dominated by the bottom squark-gluino contribution, which can enhance or suppress the Yukawa 
coupling hf, for negative or positive fi, respectively. In the following we choose At/J. < and consider 
the two possible sign combinations for At and 

To be conservative, we have chosen relatively large values for the SUSY breaking parameters, on the 
order of 1 TeV, so that some super symmetric particles may not be kinematically accessible at the LC.p^ 
However, for simultaneously large fi and Mg, the size of the A;, corrections may drive the bottom Yukawa 
coupling out of the perturbative region. Thus the gluino mass is taken as Mg = 0.5 TeV for large /i and 
Mg = 1 TeV for moderate ^. The other gaugino mass parameters are M2 = 2Mi = 200 GeV {AI2 is 
relevant for the one- loop /i — > 77 amplitude). Finally, the masses of the remaining squarks and sleptons 
are set to 1 TeV. 

We calculate the properties of the MSSM Higgs bosons for each of the benchmark scenarios using the 



program Hdecay |53|, to which we have added the A;, and A^ Yukawa vertex corrections in the two- loop 
Higgs boson squared-mass matrix |p6| and the A^ Yukawa vertex corrections in the Higgs couplings to 
bb. By comparing the partial widths, total width and branching ratios of the SM-like Higgs boson of the 
MSSM and the SM Higgs boson of the same mass, we can evaluate the sensitivity of the various observables 
for distinguishing the MSSM from the SM. We define the fractional deviations of the MSSM Higgs partial 
widths from those of the SM Higgs boson of the same mass as follows: 

j-n ITmssm - TsmI .X 
6T = , (3.1) 

1 SM 

and analogously for the branching ratios, (5BR = |BRmssm — BRsm|/BRsm- This allows us to demonstrate 
which Higgs decay quantities are the most sensitive to the non-standard nature of the Higgs boson. Later, 
we will combine these individual deviations into a variable to improve the experimental sensitivity. 



If the charged Higgs boson is Hght, then our choice of the sign of Atfi is favored by the experimental constraints on 
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b^s-f |50|,|5^1. 

^*Note that the m^-independent decoupling in the large fi and At scenario depends on ^jl/Ms and At /Ms- Lower values of 
Ms would allow for correspondingly lower values of fj, and At. 
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In the next section, we discuss the expected behavior of the Higgs BRs and partial widths in the 
MSSM, with particular emphasis on the approach to the decoupling limit In general, the couplings of 
the SM-like Higgs boson of the MSSM deviate from those of the SM Higgs boson of the same mass, except 
in the decoupling limit. 



3.1 Theoretical Expectations for Direct Higgs Couplings 

Consider first the couplings of the lightest CP-even Higgs boson h to vector bosons {V = W or Z). The 
corresponding tree- level squared-coupling normalized to the SM value is: 

4^=sin2(/?-a). (3.2) 

To a good approximation, the most important radiative correction to this result can be incorporated by 
replacing the tree-level value of a by its radiatively-corrected value [obtained in Eq. |2.3| . In the decoupling 
limit, cos(/3 — a) is given by Eq. 2.12[ It then follows that 

2 2 4-2 

JhVV ^ . _ c mz sm 

•,2 ~ ^ /|^4 



9tvv ^ ^ _ c'm% sin^ 4/? 



At large tan/3, the approach to decoupling is even faster, since sin 4/3 ~ — 4cot/3 further suppresses the 
deviation of the partial width 5r{W) = Idhww / dhsMWW ~ '^\- Contours of 6r(W) for the maximal-mixing 
and large At and /i scenarios are shown in the upper left panels of Fig. |^ and Fig. ||. The behavior in the 
no-mixing scenario is quite similar, and is therefore not shown here explicitly. 

We next consider the couplings of h to up- type fermions. At tree-level, we may write (using third- family 
notation): 

= [sin(/3 - a) + cot /3cos(/3 - q)]^ , (3.4) 



dhsMtt 

where the couplings are expressed in terms of sin(/? — q) and cos(/3 — a) in order to better illustrate the 



decoupling behavior. Based on the discussion following Eq. |2.24| , we may neglect the effects of the SUSY 
vertex corrections. Then, the dominant radiative corrections can be incorporated simply by employing 
the radiatively-corrected value of a in the tree-level coupling. Note that the approach to decoupling is 
significantly slower [by a factor of m^/m^] than in the case of the hVV coupling [Eq. |3.3| . For example. 



applying Eq. 3A to the Higgs coupling to charmed quarks, one obtains 

9hcc ^ ;l ^ cm| sin 4/3 cot /3 ^ ^^^^^ 

9hsMCC ^\ 

At large tan /3, the approach to decoupling is faster due to the additional suppression factor of cot^ /3 as 
in the case of the hVV coupling. The decoupling of 5r(c) as ruA increases is exhibited in the upper right 
panel of Figs. |l| and ||. The behavior in the no-mixing scenario is again similar and is not shown. In 
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the three benchmarks considered, the W and c couphngs to the Ught Higgs boson decouple quickly for 
increasing niA, with a somewhat slower decoupling for the c quark in the low tan/3 regime. 

Finally, we turn to the coupling of h to down-type fermions, and focus on the third-generation bb and 
T~^T~ decay modes. The approach to the decoupling limit was given for g^bf, in Eq. 2.20| [with a similar 



result for ghrr easily obtained]. For tua ^ mz (and neglecting 6hb/hb which is not tan /3-enhanced) , it 
follows that 

(3.6) 



9hbb ^ 1 _ 4cm| cos 2/3 



2 

9hsMbb '^A 



■ 2o Ah 

sm p 



1 + A, 

The approach to decoupling is again slower as compared to gnw- However, in contrast to the previous 
two cases, there is no suppression at large tan/3. In fact, since A5 oc tan/3, the approach to decoupling is 
further delayed, unless c ~ 0. Thus, we expect the greatest deviation from the SM in 5T{b) and (^r(r), and 
this is confirmed in the lower panels of Figs. || and |2[ As before, the behavior in the no-mixing scenario is 
quite similar to that of the maximal-mixing scenario and is not shown here. There is a small difference in 
the behavior of 5T{b) with respect to 5T{t) at large tan j3 due to the effect of A;,, as discussed in Sec. |^. 
Note that in the large At and /i scenario. Fig. |2| exhibits the m^-independent decoupling phenomenon, 
corresponding to the case where c ~ 0.|^ If the MSSM parameters are such that the m^i-independent 
decoupling is realized, then the experimental sensitivity to rriA is greatly compromised. The value of tan (3 
at which this decoupling occurs [Eq. 2.15(| depends slightly on the two possible sign choices for ji and At 



(for a fixed value of ^At) through the dependence of 6M.fj on ht and hh (which depend on Aj and A, 



respectively). 

So far, we have discussed the sensitivity of the Higgs couplings to the MSSM parameters. However, 
experiments at the LC will also measure the Higgs branching ratios. Since h ^ bb is the dominant decay 
mode of a Higgs boson lighter than about 135 GeV (unless the hbb coupling is anomalously suppressed), 
r(6) dominates the total Higgs width. Thus, BR(&) is not as sensitive to deviations of T(b) from its SM 
value as the BRs for other decay modes. In particular, since r(I^) quickly approaches its SM value, 
(5BR(W) ~ ^Ftot almost independently of the value of tan f3. This is illustrated in Fig. |3| for the case of 
maximal-mixing, and it is generically true in the other benchmark scenarios. 

3.2 Loop induced couplings 

The decay modes h ^ gg and /i — > 77 proceed only at the loop level. In the MSSM, SUSY particles 
and additional Higgs bosons also run in the loops. Thus the deviations of T{g) and r(7) from their SM 
values depend not only on deviations in the fermion and W pair couplings to h but also on SUSY loop 
contributions. T(g) depends mainly on the t quark loop, which has a SM-like coupling for large tua except 
^^In principle, this phenomenon could also occur in the case of no-mixing or maximal-mixing. However, in these cases, since 
5A4^2 is quite small, the value of tan f3 one obtains from Eq. 2.15 would be so large (way beyond what is plotted in Fig. ^ and 



Fig. ^) , that the mA-independent decoupling takes place in a tan P region that is no longer theoretically meaningful. 
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at small tan (3. T{-y) depends mainly on the W boson loop, which has a very rapid decoupling behavior 
with 171^. Note that, in all the scenarios considered, deviations in the t and c couplings are significant at 
low tan (3, and deviations in the b coupling are significant for all tan/3 values [see Figs. |^ and ^. 

For the loop-induced Higgs couplings, there are two separate decoupling limits of relevance. In the 
first decoupling limit, discussed often in this paper, the non-zero tree-level Higgs couplings approach their 
SM values for tua ^ The second decoupling limit applies to loop- induced Higgs couplings in the 

limit of large SUSY particle masses. In this limit, the effects of the SUSY-loops decouple, and the loop- 
induced Higgs couplings are determined by loops of SM particles. If we now take rriA ^ thz so that the 
Higgs couplings to SM particles approach their SM limits, then the resulting loop-induced Higgs coupling 
should likewise approach its SM limit. However, suppose that rriA J?> rnz but the masses of some of the 
SUSY particles that contribute to the loop-induced Higgs couplings are not too large. In this case, the 
resulting loop-induced Higgs couplings will deviate from the corresponding SM values due to the SUSY 
loop contributions, which can be a sizable fraction of the SM loop contributions. 

In the decay h — > gg, the bottom and charm quark contributions destructively interfere with the top 
quark contribution, reducing the SM amplitude by several percent. Additional deviations arise from squark 
contributions. The top squark couplings to h are 

gm? cos a grrit , . . \ ■ /o ^T\ 

where g^ arises from the so-called D-term contribution to the scalar potential.^ In addition, the minus 
(plus) sign in Eq. ^ corresponds to hiiii (/it2i2), and sin26'j = 2mjXt/(M| where M^^ < M^^. The 

top-squark contribution to the h ^ gg amplitude behaves as gj^^.^./M? for large top squark mass, where 
the 1/M~ suppression is due to the loop integral. This suppression can be partially compensated by a large 
htiU coupling. In the decoupling limit [where cos(/3 — a) ~ 0], it follows that /isina -|- At cos a = sinpXt. 
Thus if Xt is large, the third term in Eq. gives a contribution to the h ^ gg amplitude proportional to 
mtXtsm26f/AI? [i = 1,2]. Note that if \Xt\ ~ Mz and the diagonal elements of the top squark squared- 
mass matrix are of the same order, then |sin20^| ~ 1 and this contribution to the amplitude decouples 
like mt/Mf., i.e., suppressed by one power of Mj, . This should be contrasted with the case of small top- 
squark mixing angle {i.e., \ sin20^| <C 1), which arises when mfXt is small compared to the difference of 
the diagonal entries of the squark squared-mass matrix. In this case, the dominant contribution to the 
^^In principle, deviations can arise from radiative corrections due to loops of SUSY particles, but these corrections will be 

a small fraction of the corresponding SM tree-level Higgs coupling. 

^^Explicitly, go = gmz{l =F (| sin^ 6w — 1) cos2S() sin(Q + (3)/4:Cos6w- This term is independent of rrit and gives rise to 
a small coupling of h to all squark species, independent of the corresponding quark Yukawa coupling, so that the squarks of 
the first two generations also contribute to h gg. Thus, in our numerical calculations we must specify the masses of the 
squarks of the first two generations; we set Mq — Mjy = Alf^ — 1 TeV in all cases, as mentioned before. With this choice, the 
contributions of the first two generations of squarks are negligible. 
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htiti couplings comes from the second term in Eq. 3.7, and the top squark loop contribution to the h ^ gg 
amplitude decouples like mf/M? This behavior of the top squark contribution with Xt explains the 
behavior of T{g) in the no-mixing and maximal-mixing scenarios (the top panels of Fig. ^). At low m^, the 
hbb coupling is enhanced over its SM value so that T{g) is suppressed due to the destructive interference 
between the bottom quark and top quark loops. In the no-mixing scenario, the top squark contribution 
enters with the same sign as the dominant top quark loop; as rriji increases, the b quark contribution 
approaches its SM value and the top squark contribution then leads to a small enhancement of T{g) at 
large tua- Thus in the no-mixing scenario (the top left panel of Fig. ^) we see that rMSSM(9) — rsM(fl') 
passes through zero for mA_ ~ 0.6-1 TeV and reaches an asymptotic value of 1-2% in the large rriA limit. 
In the maximal-mixing scenario, the top squark loop contribution is enhanced by the large value of Xt and 
enters with the opposite sign as the dominant top quark loop; thus it leads to a further suppression of T{g) 
in addition to the suppression due to the enhanced hbb coupling at low rriA- As niA increases, then, the 
b quark contribution approaches its SM value and ST(g) reaches an asymptotic value of about 6% in the 
large niA limit, without rMSSM(5) — rsM(5) passing through zero at any value of ruA- 

In the large ^ and At scenarios, the behavior of T{g) is more complicated. At low tan /?, the behavior 
of ST{g) is dominated by the decoupling behavior of the htt coupling (see Fig. §). At large tan/3, itia- 
independent decoupling effects and Af, play an important role. For < (shown in the bottom right 
panel of Fig. 6T{g) is quite significant at large tan/3, and the deviations from the SM in this region 
are independent of niA for niA ^ 0.5 TeV. For this sign of /x, is large and negative (the two terms in 
Eq. 2.10| enter with the same sign), resulting in a significant enhancement of hf, in the large tan/? region. 
Because of this enhancement of hf,, the bottom squark loop contributions are significant in this scenario 
at large tan/3; they modify the h ^ gg amplitude by a few percent, leading to several percent deviations 
in T{g). At lower values of tua, the b quark loop also contributes due to the deviation in the hbb coupling 
(see Fig. For /x > (the bottom left panel of Fig. 5T{g) exhibits m^-independent decoupling for 
tan/3 ~ 40. For this sign of /x, the two terms contributing to Af, in Eq. 2.10| enter with opposite signs, 
leading to a partial cancellation, and the overall sign of A^ is positive. Thus A;, leads to a small suppression 
of /if,, and the bottom squark loop contribution is not significant. 

Contours of 5T{'~f) are shown in Fig. ^. In the SM, r(7) is dominated by the W boson contribution; at 
low niA ^ 0.2 TeV, the deviation of r(7) from its SM value is dominated by the deviation of the hWW 
coupling (see Figs. || and |2|). Because this one-loop decay gets contributions from all charged particles 
that couple to /i, it depends not only on the parameters of the squark and slepton sectors but also on the 



charginos. It was shown in Ref. [54| that near the decoupling regime, only top and bottom squark and 
chargino contributions can generate a sizable deviation from the SM partial width. In the no-mixing and 
maximal-mixing scenarios, we choose M2 = — /i = 200 GeV. This choice leads to large Higgsino-gaugino 
mixing, so that the couplings of h to chargino pairs are large. Thus at low tan /3 and niA ^ 0.2 TeV, the 
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chargino contribution to 6T('y) is responsible for the bulk of the deviation from the SM. As tan /3 increases, 
the Higgsino-gaugino mixing decreases and the chargino contribution becomes smaller (see the top two 
panels of Fig. ^). As in the case of /i ^ gg, there are also top squark contributions to h ^ 77; these 
are very small in the no-mixing scenario but somewhat larger in the maximal-mixing scenario due to the 
large Xt enhancement of the htiti couplings, as discussed before. In the maximal-mixing scenario, the top 
squark and chargino contributions to the amplitude enter with opposite signs. Then as tan (3 increases, 
the chargino contribution shrinks until it is the same size as the top squark contribution; this occurs for 
tan /? ~ 15 — 20, where ST{'j) goes to zero in this scenario. At larger values of tan /3, the top squark 
contribution is responsible for the bulk of the deviation from the SM. In the large /U and At scenarios (the 
bottom two panels of Fig. |^, we have M2 <^ so that the Higgsino-gaugino mixing is very small and 
the couplings of h to chargino pairs are suppressed. Thus in these scenarios the chargino contribution is 
very small, and appears only at low tan /?. For < 0, the deviation at very large tan /? and low rriA is due 
to the deviation in the hbb coupling in the b quark loop, as in the case of /i — > gg. 

4 Higgs boson measurements at future colliders 

Figs. [l|-|5| demonstrate our theoretical expectations for the behavior of partial widths and branching ratios 
in the MSSM with respect to the SM. In this section, we fold those results with the expected experimental 
resolution from the next generation of experiments at linear colliders and compare it to what will be known 
from hadron collider experiments. 

4.1 Anticipated experimental uncertainties in Higgs branching ratios and couplings 

We expect quite sensitive measurements at the LC of both the Higgs production cross sections and the 
BRs for the most important Higgs boson decay modes. Combining these measurements, it is possible to 
extract the total Higgs boson width and the partial widths for the various decay modes. The measurement 
of the total cross section for e^e^ Zh yields a measurement of g^zz- Branching ratios are determined 
by selecting Higgs events in Zh production, where Z , using the recoil mass technique and counting 

the different types of events in the invariant mass region of the Higgs boson signal. These branching ratio 
determinations are model independent, even if invisible Higgs decays are present. Finally, by including 
other Z decay modes in hZ production as well as Higgs production via vector boson fusion, one can further 
improve the statistical precision of the Higgs branching ratios and Higgs couplings. 

The expected experimental uncertainties in the measurement of BRs at the LC for a 120 GeV SM-like 
Higgs boson are summarized in Table |2|. The first row shows the results assuming 500 fb~^ of integrated 
luminosity at ^ = 350 GeV 0. At ^ = 500 GeV, the e+e" ^ Zh cross section for nih = 120 GeV 
is about a factor of two smaller than at 350 GeV ||5^; thus roughly twice as much integrated luminosity 
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Decay mode: bb WW* t^t cc gg 77 



Ref. 


2.4% 5.1% 


5.0% 


3.5% 5.5% 


19% 


Ref. Q 


2.9% 9.3% 


7.9% 


39% 18% 




Ref. ll| (scaled) 








14% 


theory uncertainty 


1.4% 2.3% 


2.3% 


23% 5.7% 


2.3% 



Table 2: Expected fractional uncertainty of BR measurements at an e^e LC for a 120 GeV SM-like Higgs 
boson. Results are shown from Ref. § (500 fb'^ at ^ = 350 GeV) (first row); Ref. [|l2| (500 fb^^ at 
^ = 500 GeV) (second row); and Ref. ^ (1 ab^^ at ^ = 500 GeV, scaled to 500 fb"!) (third row). 
The theoretical uncertainty of the predicted Standard Model branching ratios is given in the fourth row 
(see Sec. O). 



(i.e., about 1 ab^^) would be needed at 500 GeV to obtain the same statistical precision on Higgs BRs. 
To estimate the precision on Higgs BRs that can be obtained with 500 fb~^ at ^/s = 500 GeV, the results 
for = 350 GeV shown on the first line of Table |2| should be reduced by a factor of 1.5 (corresponding 
to the square root of the ratio of the corresponding Zh cross sections for mh = 120 GeV)J^ The second 
row of Table |2| shows the results of a similar study ||l^ for the branching ratios of a 120 GeV SM-like 
Higgs boson with 500 fb-i at ^ = 500 GeV.0 Finally, we consider the results of a dedicated study of 
the BR(7) measurement [13| for ^/s = 350 and 500 GeV, both without and with beam polarization (80% 



left-handed electron polarization and 40 or 60% right-handed positron polarization) chosen to enhance 
the Higgsstrahlung and WW fusion cross sections. At ^/s = 500 GeV and the highest polarizations, a 
measurement of BR(7) with an experimental uncertainty of 9.6% is possible with 1 ab^^. Scaling this to 
500 fb~^ to compare with the other studies yields a precision of about 14%, as shown in the third row of 
Table ^. Without beam polarization, this deteriorates to 16% (23%) with 1 ab^^ (500 fb^^). 

The 5r(7) deviations shown in Fig. ^ are typically too small to be observed at an e'^e~ linear collider. 
However, the operation of the LC as a photon-photon collider offers the possibility of a direct measurement 
of r(7) from the s-channel Higgs production cross section (1(77 — > h). Estimates for the precision obtainable 
on 17(77 h) X BR(/i — > bb) range from 2-10% for a light Higgs boson with mass between 120 and 160 
GeV ^from a 77 collider running at ^/s^ ~ rrih/O.S (giving peak 77 luminosity at the Higgs mass) and an 
integrated luminosity corresponding to 400 fb~^ of e^e^ luminosity [^. Combining a 2% measurement 
of cr(77 — > /i) X BR(/i bb) at a 77 collider for a 120 GeV Higgs boson with a 3% measurement of 
^^This scaling is only approximate since it does not take into account the Higgs production by vector boson fusion employed 
in the analysis of Ref. . 

^^Note the very different predictions in Table ^for the precisions of BR(c) and BK{g), which depend on very good charm 

and light quark separation. The authors of Refs. and |l2j are working together to resolve these discrepancies js^ . 
■^"The variation in precision is due to the decline of BR(/i bb) as nih increases. 
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2 2 2 2 2 2 2 

9hbb 9hWW 9hZZ 9hcc 9hTT 9hgg 9htt 



experimental uncertainty 


4.4% 2.4% 


2.4% 7.4% 6.6% 


7.4% 10% 


theory uncertainty 


3.5% 


24% 


3.9% 2.5% 



Table 3: Expected uncertainty of measurements of squared couplings (equivalently partial widths) for a 
120 GeV SM-like Higgs boson from HFitter [0,|5^, assuming 500 fb^^ at a/s = 500 GeV, except for the 
measurement of gf^^^ which assumes 1000 fb^^ at ^/s = 800 GeV. The second line shows the theoretical 
uncertainty (see Sec. |4.2| ). 



BR(/i — > bh) at the e^e~ LC, we find that r(7) could be extracted with an uncertainty of about 3.6%. As 
Fig. ^ shows, such a measurement would not exhibit a significant deviation from the SM prediction in our 
benchmark scenarios unless itia ^ 200 GeV or tan [3 is very small; however, a more detailed analysis of the 
MSSM parameter space should be performed in light of this expected precision. 

From the measurement of g\y/\/\/ based on the production cross section, the partial width r(VF) can 
be calculated and combined with the measurement of BR(VF) to determine the total Higgs width Ftot- 
The expected resolution on gf^^z about 3% yields a measurement of the total width to roughly 6% 
accuracy [@, ^] . This method for extracting the total Higgs width is more accurate than using the photon 
collider mode.^ However, the photon collider mode is still quite useful for the high precision measurement 
that it provides of the partial width r(7). In other models, such as the (non-SUSY) two Higgs doublet 
model, this measurement can be essential for distinguishing the extended model from the SM in some 



regions of parameter space [58|. 



Combining the measurements of Higgs boson BRs and production cross sections that can be obtained 
at the LC, the Higgs couplings to SM particles can be extracted [^5|. The results of a minimization 
using HFitter 0, ^ are summarized in Table ^. The first six couplings listed in Table |3| are extracted 
from h BRs into 66, WW*, cc, tt and gg and the production cross sections in the Higgsstrahlung and 
WW fusion modes, assuming 500 fb^^ at ^/s = 500 GeV [the first five of these are given in Ref. ||5^]. The 
htt coupling can be measured indirectly from the LC measurements of h ^ gg and /i ^ 77 if one assumes 
that SUSY loop contributions are negligible; however, this is a model-dependent assumption that we wish 
to avoid. A direct measurement of gf^^^ can be obtained from the e+e~ tth cross section |5^,|6^. Such a 
measurement requires running at higher ^/s = 800-1000 GeV in order to avoid kinematic suppression of 
■^^As noted earlier, r(7) can be measured to about 3-4% accuracy in 77 collisions. The total Higgs width is then extracted 
by combining this measurement with BR(7). Using the range of uncertainties in BR(7) given in Table we see that the 
uncertainty in the total Higgs width extracted by this method is dominated by the large uncertainty in BR(7). This can be 
improved somewhat by combining the LC measurements and the LHC data on r(r)/r(7) (see Table as discussed at the 
end of Sec. 4.5. Even in this case, the uncertainty in BR(7) still dominates the total Higgs width as extracted from the 77 
collider measurements. 
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the cross section; the result in Table ^ assumes 1000 fb ^ at ^/s = 800 GeV. The studies summarized in 
Tables I and I were conducted for the SM Higgs boson of mass 120 GeV, and thus are directly applicable to 
the study of a SM-like Higgs boson of the MSSM with a mass near 120 GeV, especially near the decoupling 
limit. In our benchmark scenarios, the value of rrih varies between 118 and 129 GeV. We shall assume 
that the mass dependence of the results quoted in Tables |2| and ^ is minimal in this mass range, but this 
assumption should be tested with detailed simulations. 



4.2 Theoretical uncertainties in Higgs branching ratios and couphngs 

In order to gauge the significance of an observed deviation of Higgs boson properties at the LC from the 
Standard Model expectation, one must take into account both the experimental uncertainties (statistical 



and systematic) described in Sec. 4.1 and the corresponding theoretical uncertainties for the Standard 
Model Higgs boson. Sources of theoretical uncertainty include higher order loop corrections to Higgs decay 
rates not yet computed and parametric uncertainties due to the choice of input parameters. The largest 
sources of uncertainty arise from the choice of Og, rUc and To determine the current theoretical 



uncertainty in Higgs branching ratios and couplings, we choose a<j = 0.1185 it 0.0020 [61|, mc{mc 



1.23 ± 0.09 GeV |62| and mb{mb) = 4.17 ± 0.05 GeV ]63[ (see Ref. |6J for an alternative evaluation of MS 



quark masses). By varying these input parameters in the program Hdecay |53|, we obtain the theoretical 
fractional uncertainties for the Higgs branching ratios quoted in Table For the Higgs squared-couplings 
listed in Table ^, the only significant theoretical uncertainties reside in 5^^^ and gf^^^, due to the uncertainties 
in the b and c quark masses and in (which governs the running of the quark masses from the quark 
mass to the Higgs mass). The resulting theoretical uncertainties for (7^^^ and gf^^^ (for a SM Higgs boson 
of mass 120 GeV) are 3.5% and 24%, respectively. In addition, we find a theoretical uncertainty in gf^^^ of 
3.9% due to the uncertainty in as- 

For a SM Higgs boson with rrih = 120 GeV, about 2/3 of the width is due to h ^ bb. The theoretical 
fractional uncertainties for the Higgs branching ratios to WW*, t^t^ and 77 listed in Table |2| are due 
primarily to the fractional uncertainty of the total width, which for a SM Higgs boson with mh = 120 GeV 
is mainly governed by the corresponding uncertainty in the h ^ bb width.^ The theoretical uncertainty 
in the h ^ cc decay rate is particularly large due to the relatively large uncertainty in the charmed quark 
mass. It is hoped that further theoretical work and numerical improvements in lattice calculations p5[ will 



■^■^The experimental uncertainty in (j^^^ can be reduced by combining the e'^e' tth cross section measurement with the 

measurements of h —> gg and /i ^ 77 if one assumes that SUSY loop contributions to the latter are negligible; the resulting 

uncertainty in (j^^^ is 6.0% |^ . 

^■^The observed uncertainty in mt has only a small effect on the predictions for the h —* gg and /i ^ 77 decay rates. This 

is not surprising given that the top quark mass is the most accurately known of all the quark masses! 

■^''For larger values of the Higgs mass, the /i — > 66 branching ratio is smaller and the uncertainty in the total width, which 



is now dominated hy h ^ WW^*\ is correspondingly reduced. 
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reduce this theoretical uncertainty by the time that LC data is available. 

Finally, a scan of the tt threshold at the LC will yield a value of nit with an uncertainty of about 
100 MeV Thus, the theoretical error expected for the Standard Model Higgs coupling to ti due to the 
top-quark mass uncertainty will be negligible. The remaining uncertainty in gf^^^ is due to uncalculated 
higher order QCD corrections to the e~^e~ — > tth cross section. We estimate this uncertainty to be about 
2.5% based on the renormalization scale dependence in the next-to-leading order (NLO) QCD result for 
nih = 120 GeV and = 1 TeV |67]. 



4.3 Branching ratio analysis 

A number of the Higgs BRs can be measured to higher accuracy than the total Higgs width. Thus BR 
measurements alone are valuable for distinguishing the SM Higgs boson from a MSSM Higgs boson. To 
illustrate the potential of the LC, contours of (Eq. |3.1| ) are shown in Fig. ^ over the niA — tan/3 
plane for the benchmark scenarios. Contours of (^BR(5) = 3 and 6% and (5BR(14^), 6mi{g) = 8 and 
16% were chosen, corresponding roughly to one and two times the expected experimental uncertainties 
quoted in Table or approximately one and two sigma deviations from the SM.^ As shown in Table 
the theoretical uncertainties in these three BRs for the SM Higgs boson are smaller than the expected 
experimental uncertainties. 

In the four panels of Fig. ^, the solid, long-dashed, and short-dashed lines are contours of (5BR(6), 
(5BR(1^) and 5BIl{g), respectively. Although 5T{b) is quite large over much of the parameter space, 
(5BR(6) is smaller because the increase in T(b) also significantly increases Ftot- Because (5r(l^) quickly 
approaches zero for increasing mA_, 6mi{W) indicates variation in the total Higgs width, and is more 
sensitive than 5BR(6), except for the case of maximal mixing. In regions of parameter space where ST(g) 
approaches zero (see Fig. 63K{g), like JBR(H^), is sensitive to variations in the total width. 

For the maximal-mixing scenario, the mass of the SM-like Higgs boson near the decoupling limit is 
roughly 10 GeV heavier than in the other benchmarks (see Table |l]), so that the relative contribution 
of r(6) to Ftot is decreased. Therefore, deviations in ghbf, are not as diluted in the BR measurement as 
in the other scenarios, and the measurement of (5BR(6) yields superior sensitivity at large tan /3, around 
rriyi <^ 600-700 GeV at 2a. One should interpret this result with caution, however, since the accuracies 
for BR measurements are based on the simulation of a 120 GeV SM Higgs boson. In the maximal-mixing 
scenario, BTi{g) deviates by more than 8% from its SM value for niA ^1-4 TeV. At 2a the reach in 6JiIi{g) 
is roughly rriA ^ 600 GeV. In the no-mixing scenario, 6BK{g) and (5BR(H^) give comparable reach in niA', 
at 2a the reach is rriA ^ 425 GeV. For comparison, in the no-mixing scenario deviations in BR(6) yield 
sensitivity at 2a for niA ^ 300 GeV for tan/? > 5. 



■^^We choose typical values rather than using the precisions quoted in the experimental studies |7|,|l2| because we expect the 
exact results of these studies to change as the LC detector design evolves and experimental techniques improve. 
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The large /x and At scenario demonstrates the complementarity of the LC and the hadron colliders. 
For tan/? < 20, where the heavy MSSM Higgs bosons can be missed at the LHC, Jili{g) gives the greatest 
reach in m^, allowing one to distinguish the MSSM from the SM Higgs boson at 2a for tjia ^ 350-450 
GeV, depending on the value of tan p. At larger values of tan /?, the large fi and At scenarios have regions 
of m^i-independent decoupling where the SM-like MSSM Higgs boson cannot be distinguished from the 
SM Higgs boson even for very low values of rriA- In fact, in these scenarios it is possible for h to be 
indistinguishable from the SM Higgs boson at the LC, while at the same time niA < 250 GeV so that the 
heavy Higgs bosons will be directly observed at a 500 GeV LC through e^e~ HA, H^H~ . Moreover, 
in our scenarios with Atfi < 0, for ^ > < 0] the m/^-independent decoupling occurs for tan /? ~ 40 
[tan/? ~ 33] .0 For such large values of tan/?, the heavy MSSM Higgs bosons would be discovered at the 
LHC even for niA above 500 GeV Note also that for large fj, and At and large tan/3, the correction 
A;, is quite large, and modifies the b quark Yukawa coupling from its SM value. The effect of A;, on the 
couplings of the heavy MSSM Higgs bosons does not decouple for uia ^ mz (i.e., for tan a tan/? = — 1). 
Thus A5 could have a significant effect on the discovery of the heavy Higgs bosons at the Tevatron and 
LHC in this region of parameter space by modifying their production cross sections and decay branching 
ratios. In particular, the value of tan f3 extracted from the A and H production rates at the LHC could 
not be unambiguously determined without knowledge of the value of A^. LC data would then be of great 
value for disentangling the Af, dependence (see Sec. ^) so that the value of tan /? extracted from heavy 
Higgs boson measurements can be compared to the value obtained from other sectors of the theory. 

Clearly, from Fig. ^, the regions of the itla — tan P plane in which the MSSM and SM Higgs bosons can 
be distinguished from one another depend strongly on the supersymmetric parameters, and the sensitivity 
comes from different measurements for different sets of MSSM parameters. In isolation, measurements 
of one or two sigma deviations in individual BRs would not be a significant probe of the MSSM. In 
combination, however, the measurements are much more powerful, as indicated by the analysis presented 
in the next section. 



4.4 analysis of couplings 

In order to make a quantitative assessment of the ability of the LC to discriminate between the SM-like 
Higgs boson of the MSSM and the SM Higgs boson, we combine several observations and compute the 
compatibility with the SM using a test. In particular, = Si(^i^^^^ ~ where z is a 

decay product and X is a BR, F, a Higgs squared-coupling or a ratio of these quantities. The values 
include the experimental resolution and any theoretical uncertainty on the quantity Xi. Motivated by 



the HFiTTER results [55|, we choose the Xi to be the squared couplings of Higgs bosons to various final 
■^^As explained in Sec. 3.1, the value of tan/3 at which this mA-independent decoupling occurs has a small dependence on 
the signs of fi and At- 
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states. The significance of a particular value for ^ depends on the number of observables that have been 
combined. Our results indicate that the addition of several variables does not necessarily improve the 
significance. For example, the /177 coupling is not measured very well, and the relative theoretical error 
in the predicted hcc coupling is too large. Moreover, the hWW and hcc couplings quickly approach their 
SM values for increasing niA- Thus measurements of these three couplings do not add to the significance 
of our results. Therefore, we compute the combining bb, t~^t~' and gg squared-coupling measurements, 
adding in quadrature the experimental and theoretical uncertainties in each squared-coupling [see Table . 

The results are shown in Fig. ^ for the benchmark scenarios with contours corresponding to 
68,90,95,98 and 99% confidence levels. For tan/3 ^ 5, the no-mixing and maximal-mixing scenarios 
can be distinguished from the SM at 95% [99%] confidence level for tua ^ 600 [500] GeV and 650 [600] 
GeV, respectively. The large At and fi scenarios have regions at some large tan /? values that are indistin- 
guishable from the SM for any value of rriA- For tan/3 < 20, however, the large fi and At scenarios can 
be distinguished from the SM at 95% [99%] confidence level for niA ^ 450 [425] GeV. Note also that for 
/i < 0, the large A;, effects lead to significant deviations from the SM even for ^ 2 TeV at the largest 
values of tan /?. 



4.5 Complementarity to hadron collider measurements 

Measurements of Higgs boson properties will be available from the Tevatron and the LHC. This data will 
most likely be available before the LC is operational. Based on our current understanding of the MSSM 
Higgs boson properties and experimental capabilities, these hadron colliders will observe a light, SM-like 
Higgs boson and, perhaps, other non-SM-like Higgs bosons if their couplings to heavy flavor are enhanced 
over the SM {e.g., H and A will be observed in bbH/A if tan/3 is sufficiently large [@,^). 

If a light, SM-like Higgs boson is discovered, various combinations of production cross sections times 
branching ratios will be measured at the LHC to about 10-20% ||6^ assuming an integrated luminosity of 
100 fb~^ at each of the two detectors. The uncertainties on most of these measurements are dominated by 
statistical error. From these measurements, various ratios of the partial widths of a 120 GeV SM-like Higgs 
boson to ZZ*, WW*, 77, t~^t~ , and gg can be extracted with uncertainties between 15 and 30%. Expected 
uncertainties from Ref. [S8| are summarized in Table |^,^ along with the corresponding uncertainties from 
LC measurements with 500 fb~^ at ^ = 350 GeV [0] and from HFitter at the LC with 500 fb"^ 

at ^ = 500 GeV. 

Assuming that the hWW and hZZ couplings are related by the usual SU(2) relation, that r(6) is related 
to T{t) by the SM relation, and that only the decay modes bb, ZZ* , WW*, 77, t^t~ , and gg are needed 
to estimate the total Higgs width, then r(VF) and Ftot can also be extracted from LHC measurements 
■^^ The decay modes to 77, ZZ* and WW* were considered for inclusive Higgs production (dominated by gluon fusion) and 
the decay modes 77, t'^t^ and WW* were considered for Higgs production through vector boson fusion. 
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nz)/T{w) n^)/T{w) nT)/T{w) r(r)/r(7) ng)/T{w) r(6)/r(Ty) 



LHC PI 



16% 



15% 



15% 



15% 



LC 



HFlTTER [55 1 



20% 



3.4% 



7.1% 
6.6% 



20% 



71 



5.6% 
2.4% 



Table 4: Expected uncertainties for the ratios of partial widths of a 120 GeV SM-like Higgs boson from the 
LHC II (100 fb"^ at each of the two detectors), the e+e" LC §] (ratios of BRs; 500 fb^^ at ^ = 350 
GeV), and the e+e" LC using HFitter 0,|5|,||] (500 fb^^ at ^ = 500 GeV). 

with uncertainties of 10 and 20%, respectively |^8|. (The second assumption is violated by the Af, Yukawa 
vertex corrections; this affects the determination of both r(VF) and Ftot-) The branching ratio for the 
decay to bh is more difficult to measure at the LHC due to QCD backgrounds. T{h) can be extracted to 
about 50% from vector boson fusion Whjj events [^. 

From Table ||, it is clear that LC measurements would give a significant improvement over the hadron 
collider measurements. The LC also offers the advantage of model-independent measurements of the 
Higgs boson branching ratios, even if invisible Higgs boson decays are present, and a model-independent 
determination of the Higgs total width. Generally speaking, except for the rare 77 decay, the LC has 
measurement uncertainties that are smaller by at least a factor of 2 compared to the LHC. The precision 
of bb measurements are improved by an order of magnitude at the LC. 

Some information may also be available from the hadron colliders regarding niA or m fj± and tan /3 Q . 
For example, if the SUSY spectrum is such that Af, and A^ are negligible, then tan (3 can be measured 
at the LHC from bbH/A, H/A — > r"*"T~ event rates with a 5-25% statistical uncertainty and about a 
10% luminosity uncertainty |^,^. However, as indicated earlier, the extraction of tan/3 from LHC Higgs 
boson measurements will require knowledge of A^ and A,-. Moreover, several interpretations of high- mass 
excesses of r's may exist, especially within the MSSM at large tan /3. As we have discussed already, in 
most regions of parameter space the LC can provide an indirect measurement of niA for values of tjia 
significantly beyond the e~^e~ AH kinematic limit. 

From Tables ^ and ^ we see that the combination of LHC and LC measurements will significantly 
improve our knowledge of BR(7) in the absence of a 77 collider. In particular, the ratio of Higgs par- 
tial widths F(r)/F(7) can be determined to 15% accuracy using ratios of production cross sections times 
branching ratios in Higgs production through vector boson fusion. Combining this with the LC measure- 
ment of BR(t) to 5-8%, we find that BR(7) can be extracted with a precision of 16-17%. If the BR(7) 
measurement at the LC of 14-19% is combined with this BR(7) extraction from LHC data, a "world 
average" precision of 11-13% can be obtained. 
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5 Extracting SUSY parameters: A;, 



In this section we examine the possibihty of extracting MSSM parameters from measurements of the 
properties of the SM-hke Higgs boson. We concentrate on the SUSY vertex corrections to the hbb couphng, 
Ah (see Sec. |2|). As noted in Sec. |4.3| , when A;, is large it can have a significant effect on the phenomenology 
of the heavy MSSM Higgs bosons. In particular, at large tan /3, knowledge of Ah will be needed to extract 
the value of tan (3 from measurements of the A and H production rates at the LHC Because Ah depends 
on a combination of MSSM parameters (Eq. |2.10| ), a measurement of this quantity provides a constraint on 
MSSM parameter space. However, additional information is needed to fully disentangle these parameters; 
e.g., the gluino and bottom squark masses can be obtained from measurements at hadron colliders. To 
extract Ah directly from the data in the large tan /3 regime, we shall consider the following ratio of Higgs 



couplings [see Eq. |2.26 | 

9hbb — 9hTT 9hhb — Qhrr 



-Afe, (5.1) 



Qhtt — Qhbb ghcc — Qhbb 

where we have employed the notation introduced at the end of Sec. ^ for Higgs-fermion couplings normalized 
to their Standard Model values. 

At the LC the hhh and hrr couplings can be extracted from measurements of BR(6), BR(r) and the 
total width of /i, as discussed in Sec. In the up-type quark sector, the measurements of h couplings to 
charm- and top-quark pairs are possible at the LC. In h decays, we have access to the hcc coupling through 
BR(c). This measurement, like those of the hhb and hrr couplings, can be made in the initial stage of 



the LC running at ^/s = 350 or 500 GeV. However, as discussed in Sec. 4.2, the theoretical uncertainty in 
QhsMcc is very large due to the large uncertainty in the charm quark mass. If the uncertainty in the charm 
quark mass is not significantly reduced by the time that LC data is available, then the resulting theoretical 
uncertainty in g^cc would dominate the uncertainty in the value of A;, extracted from Eq. |5.1| . In contrast, 
the theoretical uncertainty in the htt coupling is very small due to the precision in the top quark mass 



measurement. As discussed in Sec. 4.1, the coupling ghu can be measured directly from the e+e" — > tth 
cross section; however, this requires LC running at a higher energy, e.g., ^/s = 800 GeV. 

The LC measurements of h decay properties are sensitive to the squares of the Higgs couplings ^^^j 
through the h partial widths; therefore the signs of the normalized couplings ghff in Eq. are not 
measured directly in the corresponding BRs. This leads to a potential four-fold ambiguity in the extraction 
of \ Ah\. However, LC measurements oi h ^ gg and h ^ are sensitive to the relative signs of the Higgs 
couplings to the particles running in the loops; for example, flipping the sign of ghtt while leaving that 
of 9hww fixed leads to a large deviation in the h'j'j coupling from its SM value, while the hgg coupling 
is sensitive to the relative sign of ghtt and g^hb (assuming that SUSY loops are not important). Near the 
decoupling limit, the radiative corrections to the CP-even Higgs mixing angle a are not so large as to 
change the sign of tana from its tree- level negative value, and |At-| < \Ah\ < 1. In this case the signs of 
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Qhff are the same as in the SM, so all the ghff in Eq. 5.1 are positive. 

In the discussion above, we have assumed that h is SM-like, which is appropriate once niA is sufficiently 
above m™'^-^ (in practice, rriA ^ 150 GeV is usually sufficient). For values of ttia ^ m™^^ and tan/? ^ 1, 
it is H which has SM-like couplings to vector boson pairs (i.e., the roles of h and H are interchanged). 
However, a separate analysis is not required since we also have [see Eq. pel 

9Hbb — 9Htt ^ QHbb ~ gHrr ^ 2^ 

gmt — gmb gucc — gmh 
The extraction of the normalized couplings gnff follows the same procedure as outlined above. 

The upper two panels of Fig. ^ show the fractional error in the determination of from measurements 
of the Higgs couplings to 66, r+r^ and tt, for the benchmark large and At scenario (see Table Q), in 
which Afe is quite sizable. The error in Af, is calculated using the fractional uncertainties in the Higgs 
couplings given in Table ^ In the two upper panels of Fig. ^ in the region outside of the dot-dashed 
contours [and also within the smaller closed dot-dashed circles], the error on Af, is larger than A;, itself, 
and no meaningful value for A^ can be obtained. Three other contour lines are shown corresponding to a 
fractional uncertainty in A;, of 50% [dotted], 10% [dashed] and 5% [solid]. Within the region bracketed only 
by solid contours, the relative error on A;, is less than 5%. There is an approximate (distorted reflection) 
symmetry in these contours about ruA = m^^^ = 119 GeV. This corresponds to the interchange of roles 
between h and H at large tan [5 as noted above. Note also the convergence of the contours at tan /3 ~ 41 
(left panels) and tan /3 ~ 33 (right panels) corresponding to the tan (3 values at which the myi-independent 
decoupling occurs. 

For positive /x (the upper left panel of Fig. ^), A;, can only be distinguished from zero at the 2a level 
(i.e., a 50% measurement) for moderate to large tan [3 and tua ^ 160 GeV. For such low values of m^, the 
heavy Higgs bosons will be directly observed at the LC. In contrast, for negative ^ (the upper right panel 
of Fig. ^; note the change in the horizontal scale), Af, can be determined with 10% accuracy even for tua 
as large as 500 GeV for large tan/3. This significant difference in sensitivity to A^ between the positive 
and negative cases is a combination of two effects. First, for negative the bottom-squark-gluino and 
top-squark-Higgsino contributions to Af, enter with the same (negative) sign, while for positive /x they 
enter with opposite signs (see Eq. p.lOj ). Thus for fixed values of and |/i| (with At^ < 0), |Af,| is 
larger for negative /x than for positive fi. Second, Af, enters the 6 quark Yukawa coupling non-linearly as 
1/(1 + A;,), so that Af, < has a more pronounced effect. 

The upper right panels of Figs. || and |2| demonstrate that r(c) oc gf^^^ is very close to its SM value once 
ruA ^ 150 GeV and tan/3 ^ 5. In this region of parameter space, | cos(/? — a)| ^ 1 and 



ghcc — ghtt 



1 + cot /3 cos(/3 -a) + (cos'^{(3 - a)) ~ 1 . (5.3) 



Thus a value of Af, can be extracted by taking ghcc = ghtt = 1 in Eq. |5.1| . To the extent that this is a good 
approximation, there is no sensitivity in the value of A^ to the theoretical uncertainty in the charm or top 
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quark Yukawa couplings. Likewise, for niA ^ m™^^ and tan/? ^ 1, sin(/3 — a) is close to zero so that 

gHcc = 9Htt 1 - cot /? sin(/3 - a) + O U\r? (/3 - a)) ~ 1 , (5.4) 



and At can be extracted by taking gucc = 5/^** = 1 in Eq. 5.2. Results are shown in the lower two panels 
of Fig. ^ With the assumption that the Higgs coupling to ti is given by its SM value, the precision in 
is improved somewhat, especially for large values of tan/3 and rriA and fi < [the lower right panel of 
Fig. M. Note that in the region of the iriA — tan/3 plane in which cos(/3 — a) takes on some intermediate 



value between zero and one, neither Eq. 5.3 nor Eq. is a good approximation, and the results in the 
lower two panels of Fig. ^ are not reliable. 



6 Conclusions 

There are significant regions of MSSM parameter space where the Tevatron or LHC can discover a light 
Higgs boson which is indistinguishable from the SM Higgs boson using only hadron collider measurements. 
Additional Higgs bosons may be observed, but further measurements will be necessary to constrain the 
underlying Higgs sector. Also, there is a region of MSSM parameter space for moderate tan /3 and moderate 
to large rriA where only a single light Higgs boson with SM-like couplings would be observed in experiments 
at the Tevatron, the LHC, and a y/s = 350 or 500 GeV e~^e~ linear collider. In this case, precision 
measurements of the properties of the light Higgs boson at the LC can be used to distinguish between the 
SM and MSSM, and further to extract or constrain the model parameters. At the LC, the production cross 
sections and the most important branching ratios of a light SM-like Higgs boson can be measured with 
precisions between 2 and 15%. We have examined how these measurements can be used to distinguish an 
MSSM Higgs boson from the SM Higgs boson. We have also demonstrated the extraction of the important 
SUSY quantity A^ that parameterizes the SUSY radiative corrections to the relation between the b quark 
mass and its Yukawa coupling. 

In general, the couplings of the lightest CP-even MSSM Higgs boson h to pairs of SM particles approach 
their SM values in the decoupling limit of large niA- Thus, the deviation of h partial widths from their SM 
values goes to zero in the limit of large m^, except for the loop-induced couplings of h to photon or gluon 
pairs, which approach their SM values in the limit that both and the SUSY masses are taken large 
simultaneously. As one moves away from the decoupling limit, the radiative corrections to the CP-even 
Higgs mixing angle a, the Yukawa vertex corrections A5, and contributions of SUSY particles to the loop 
induced hgg and hjj couplings, all contribute to deviations of the h couplings from the corresponding SM 
values. The sensitivity of these couplings in the — tan/3 plane for distinguishing the SM-like MSSM 
Higgs boson from the SM Higgs boson depends strongly on the MSSM parameters. In order to illustrate this 
strong dependence, we have chosen three benchmark MSSM scenarios that lead to very different behavior 
of the Higgs couplings. 
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The decoupling behavior of the h partial widths at large and moderate to large tan P is summarized 
qualitatively as follows: (5r(I^) ~ cot^ Sr{c) ~ m| cot^ and 6T{b) ~ ST{t) ~ m'^/m\. 
Moreover, there is a contribution to 6T{b) proportional to A;, [see Eq. which is enhanced at large tan /3 
and can further slow the approach to decoupling. Thus we expect that r(VF) will decouple very quickly 
with increasing rriA (so that 5BR(VF) ~ STtot over most of the parameter space), r(c) will decouple more 
slowly at low tan /5 but quickly at large tan (3, and r(6) and r(r) will decouple quite slowly, at all values of 
tan /3. This expected behavior is borne out in the no-mixing and maximal-mixing scenarios. In the large ^ 
and At scenarios the decoupling behavior of these partial widths is more complicated, because the radiative 
corrections to the CP-even Higgs mixing angle a are quite significant in some regions of the niA — tan /3 
plane. For some (large) values of tan/3, cos(/3 — a) is driven to zero by the radiative corrections to a, 
leading to "m^i-independent decoupling". In particular, r(6) and r(r) attain their SM values at much 
lower values of than would be expected from the tree-level decoupling behavior. 

The myi-independent decoupling is not an exclusive feature of the specific large values of /i and At 
chosen in our benchmark scenarios. Additional sets of MSSM parameters that exhibit m^-independent 
decoupling can be found by scaling fi, At, Mg and Ms by a common factor; in this case the value of tan f3 
at which the myi-independent decoupling occurs varies primarily due to the log(M|/m^) dependence of 



[see Eq. 2A\. In particular, we have found sets of MSSM parameters with relatively small Ms that 
exhibit myi-independent decoupling even for tan f3 less than 10, within the parameter region in which the 
heavy MSSM Higgs bosons would be missed at the LHC. In this case, due to the relatively low value of 
Ms, deviations from the SM Higgs couplings may be observable in /i — > or /i ^ 77 due to top-squark 
and bottom-squark loop contributions. 

In the case of the loop-induced Higgs couplings to gluon or photon pairs, deviations from the SM 
widths can occur due to SUSY loop contributions even in the limit of very large niA. These SUSY loops 
decouple at large SUSY particle masses, but the decoupling can occur more slowly if Xt is large (as in the 
maximal-mixing scenario), resulting in several percent deviations of T(g) from its SM value even for TeV 
mass squarks. 

Because T{b) dominates the total width of a 120 GeV SM-like Higgs boson over most of the MSSM 
parameter space, a change in T{b) causes a corresponding change in Ttot, so that the resulting change in 
BR(6) is relatively small. In contrast, branching ratios for rarer Higgs decays such as BR(l^) and BR(g') 
are sensitive both to changes in r(VF) and T(g), respectively, and to changes in Ttot- Thus, although BR(?)) 
will be the best-measured BR and that r(6) typically deviates significantly from its SM value, BR(VF) and 
BR((;) will in general be more sensitive to deviations from the SM than BR(6), making their high-precision 
measurements a priority. 

To demonstrate how LC measurements of Higgs properties can be used to distinguish a SM-like MSSM 
Higgs boson from the SM Higgs boson, we studied the expected uncertainties in the branching ratios 
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of a 120 GeV SM Higgs boson to bb, WW* and gg, and performed a analysis using the anticipated 
uncertainties in the squared couphngs of the SM Higgs boson to bb, t'^t" and gg. In particular, we draw 
the following conclusions for our three benchmark scenarios from the analysis (a somewhat lower reach 
in niA is obtained using the analysis of individual BR measurements) : 

• In the no-mixing and maximal-mixing scenarios, the MSSM Higgs boson can be distinguished from 
the SM Higgs boson at the 95% confidence level for niA ^ 600 and 650 GeV, respectively, for 
tan /? > 5. These limits become slightly lower for values of tan /? < 5. 

• In the large fi and At scenario, there are regions of parameter space at large tan/? in which the 
MSSM Higgs boson cannot be distinguished from the SM Higgs boson using BR and cross section 
measurements for any value of niA. For tan /3 < 20, the MSSM Higgs boson can be distinguished 
from the SM Higgs boson at the 95% confidence level for niA ^ 450 GeV. 

The precise reach in niA in the maximal-mixing scenario should be viewed liberally, since the SM-like Higgs 
boson mass can reach values near 130 GeV, and some of the coupling and BR measurements used in this 
study have poorer precision for ruh > 120 GeV. 

Finally, we have shown that when A(, is sufficiently large, information on can be obtained from 
the measurements of the couplings gh^b, g^rr aiid ghtt [oi' S/icc] unless the MSSM parameters are such that 
m^i-independent decoupling occurs. The precision achievable in this measurement depends on the size and 
sign of A5; these depend in part on the signs of fiMy and j^iAt- Combining precision measurements in the 
Higgs sector, such as those leading to the determination of A5, with direct measurements of the SUSY 
spectrum will be of great value in deciphering the underlying super symmetric structure of the theory. 
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Figure 1: Deviations of Higgs partial widths from their SM values in the maximal-mixing scenario. 
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Figure 2: Deviations of Higgs partial widths from their SM values in the large /x and At scenario, with 
At = -^i = 1.2 TeV. 
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Figure 3: Deviations of BR(VF) and Ftot from their SM values in the maximal-mixing scenario. 
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Figure 4: Deviations of the partial width T{g) from its SM value in the no-mixing scenario (top left), the 
maximal-mixing scenario (top right), and the large ^ and At scenario with fi = —At = 1.2 TeV (bottom 
left) and jj, = —At = —1.2 TeV (bottom right). In the no-mixing scenario, r(gi)MSSM — r(fi')sM changes 
sign and thus passes through zero [i.e., 6T{g)=0] along a contour between the two ST{g) = 0.01 contours 
(solid lines) exhibited in the top left panel. 
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Figure 5: Deviations of the partial width r(7) from its SM value in the three benchmark scenarios. 
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and 16% (short-dashed) in the three benchmark scenarios. 
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Figure 7: Contours of for Higgs boson decay observables in the benchmark scenarios. The contours 
correspond to 68, 90, 95, 98 and 99% confidence levels (right to left) for the three observables g\i,^, 9hTTi 
and gl . 
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Figure 8: Contours of the fractional uncertainty in the determination of Af, in the large // and At scenario 
with /i = -At = 1.2 TeV (A5 > 0) (left) and /i = -At = -1.2 TeV (A;, < 0) (right). The upper two 
plots use the Higgs couplings to 66, t~^t~ and ti as inputs while the lower two plots use only the first two 
couplings as input while taking the Higgs coupling to ti equal to its SM value (see Eqs. |5.3| and 5.4). In 
the = -At = 1.2 TeV scenario (left), we plot 0.1 TeV < niA < 0.25 TeV, while in the = -At = -1.2 
TeV scenario (right), we plot 0.1 TeV < rriA < 0.75 TeV. 
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